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ABSTRACT 

A new 7-beam 6—7 GHz receiver has been built to survey the Galaxy and the Magel- 
lanic Clouds for newly forming high-mass stars that are pinpointed by strong methanol 
maser emission at 6668 MHz. The receiver was jointly constructed by Jodrell Bank 
Observatory (JBO) and the Australia Telescope National Facility (ATNF) and allows 
simultaneous coverage at 6668 and 6035 MHz. It was successfully commissioned at 
Parkes in January 2006 and is now being used to conduct the Parkes- Jodrell multi- 
r% \ beam maser survey of the Milky Way. This will be the first systematic survey of the 

■ entire Galactic plane for masers of not only 6668-MHz methanol, but also 6035-MHz 

excited-state hydroxyl. The survey is two orders of magnitude faster than most pre- 
vious systematic surveys and has an rms noise level of ~0.17 Jy. This paper describes 
the observational strategy, techniques and reduction procedures of the Galactic and 
Magellanic Cloud surveys, together with deeper, pointed, follow-up observations and 
complementary observations with other instruments. It also includes an estimate of 
the survey detection efficiency. The 111 days of observations with the Parkes telescope 
have so far yielded >800 methanol sources, of which ~350 are new discoveries. The 
whole project will provide the first comprehensive Galaxy-wide catalogue of 6668-MHz 
and 6035-MHz masers. 

Key words: stars: formation, Masers, Surveys, (galaxies:) Magellanic Clouds 



1 INTRODUCTION 

The theory and understanding of the formation of high-mass 
stars is a fundamental unsolved problem in astronomy. Cur- 
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rent star formation models encounter diffi culties in creating 
stars much more massive than ~8 Mq fe.g. lKetoll2003l ; lYorkel 
|2004) and observations of these events are hampered by the 
speed with which the process occurs, and hence the relative 
rarity of objects in the act of accreting t he majority of their 
final mass. However, since its discovery l|Mentenlll99ll ). the 
6668-MHz methanol maser has been recognised as one of the 
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clearest signposts to the formation of high-mass young stars. 
It is the second strongest cosmic maser known, surpassed 
only by water (H 2 0) at 22 GHz. The 6668-MHz methanol 
maser is widespread and, unlike all other strong masers (OH, 
H2O and SiO), is only f ound close to high- mass young stars 
(e.g. iMinier et al.ll2003h and sometimes associated with ul- 
tra compact (UC) Hn regions. Additionally, the methanol 
mase r s may trace earlier hot core phases (e.g. IMinier et al.l 
120(4] 

Purc ell et al.ll2006|). al though recent kinematic work 
by Ivan der Walt et al.l 1 20071) argues otherwise. It is possi- 



ble that heating of the hot molecular core, or shock heating 
from an outflow liberates methanol from dust mantles and 
so provides t he high column density needed to give strong 
maser action (|Cragg et al.ll2002l ; [Codella et al.ll2004 ). 

During the past decade, extensive 6668-MHz methanol 
maser searches have been undertaken in our Galaxy using 
two different strategies: (1) targeted searches toward colour- 
selected infrared (IR) sources and known regions of intense 
star formation (e.g. OH and H2O masers); (2) unbiased 
surveys covering portions of the Galactic plane in both the 
northern and southern hemispheres. These methods have 
yielded different results in terms of detection stati stics. 
Except for the first targeted survey b ylMentenl (1 199lf) and 
the survey of OH maser positions by ICaswell et al.l (119951 ) 
(with success rates exceeding 80 per cent, and further 
augmented by the discovery of offset and clustered maser 



duccd detection rates of about 15 per cent (MacLeod et al. 


1992 


: iMacLeod & Gavlardl 1 19921: iGavlard & MacLeod 


1993 


: Schutte et al.l 19931: van der Walt et al.1 1995, 


1996; 


Walsh et al.| 1997: Slvsh et al.l 19991: Szvmczak et al.1 


2000; 


Ellingscn 2007]). Generally, the targeted surveys 


were 



restricted to known or suspected sites of star formation 
and missed masers arising in regions where they were 
not expected. The unbiased searches have produced 
detection rates lying between 0.22 sources degree -2 
|Pestalozzi. Minier fc Boothl 120051 ) and 



more than 1 

sources degree"' " ifcaswelll Il996bl lal: lEllingsen et~ai1 1 19961 ). 

A total of 140 degree 2 have been surveyed to date. This 
includes a region of 16 degree 2 in the first Galactic quadrant 
which have been co vered in several unbi ased surveys (as 
evident in Table 3 of iPestalozzi et al.|[2005l ). More than 500 
sources are listed in the compilatio n of 6668-MHz methanol 
sources by IPestaloz zi et al.l l|2005l ). Fig.[T] shows how these 
masers are distributed in the Galaxy (for | b |^ 2.0°). 

Almost 50 per cent of the detected obje cts of the untar- 
;eted surveys of lEllingsen et al.l (|l996l ) and lSzvmczak etafl 
20021 ) had no previously known counterpart, clearly show- 
ing the limitation of targeted observations. In both sam- 
ples, possible counterparts to the masers have a wide 
range of infrared colours, with many conspicuously falling 
outside the UC H ll regi on IRAS selection criteria of 
IWood fc Churchwelll (|l989l ). Therefore, any survey based 
solely upon long-wavelength infrared colour selection criteria 
is likely to underestimate the number of masers, and searches 
using targets based on mid-infrared colour-selected sources 
of the Spitzer Galactic Legacy Infrared Mi d-Plane Survey 
Extraordinaire (GLIMPSE) infrared survey (|Ellingsenll20071 ) 
also seem to be inefficient. 

On the other hand, the existing unbiased surveys are 
deficient because of their incomplete Galactic coverage and 
heterogeneous sensitivity limits. A recent analysis of the ex- 



isting surveys suggests that the actual number of methanol 
masers sources in our Galaxy is about twice that curre ntly 
known (Ivan der Walt 2005; Pa ndian fc Goldsmith! 120071 ). A 
potential doubling of the number of known young high mass 
stars clearly has important implications for our view of high- 
mass star formation in our Galaxy, and provides the oppor- 
tunity of enhancing our understanding of the nature and 
origin of 6668-MHz methanol masers. 

The only way to determine the actual number of 
methanol maser sources in the Galaxy, and identify the high- 
mass star formation regions which they trace, is to carry 
out a uniform, unbiased survey of the whole Galactic plane. 
With sufficient sensitivity, such a survey will detect the vast 
majority of such sources. Since the high brightness of masers 
allows their positions and velocities to be determined to 
high precision, such a Galaxy- wide survey will provide an 
excellent probe of the global distribution of star formation 
throughout the Galaxy. In addition, the unbiased, homoge- 
nous catalogue that such a survey can produce will provide 
a statistically complete sample for future studies of the early 
stages of high-mass star formation. 



2 THE MMB GALACTIC PLANE SURVEY 

The need for a new survey of the Galactic Plane for the 
6668-MHz methanol masers led to the development of the 
Methanol Multibeam (MMB) survey. The survey goal is 
to produce a high-reliability unbiased survey of 6668-MHz 
methanol masers in the Galactic plane, with a uniform sen- 
sitivity completeness. It will allow a complete census of the 
high-mass star formation taking place in our Galaxy. In this 
way, and in combination with other Galaxy-wide surveys, 
it will be possible to accumulate the statistics necessary to 
determine the relative lifetimes of the various stages of high- 
mass star formation, and hence be able to build an empirical 
model of the evolutionary process. This paper and its fore- 
runner (I Green e t al. 2008) represent the first two in a series 
on the MMB survey, soon to be followed by the first Galactic 
results. 



2.1 Survey Parameters 

To cover the whole Galactic plane in longitude the MMB 
team has used the Parkes 64-m radio telescope in the south- 
ern hemisphere and currently plans to use the Lovell 76-m 
radio telescope at Jodrell Bank in the northern hemisphere. 
The survey area has been split, with longitudes —174° < I < 
60° observed from Parkes and 60° < I < 186° planned for 
the northern hemisphere observations. 

The latitude dist ribution of methanol masers in the 
IPestalozzi et al.1 j|2005l ) compilation has a FWHM of ~0.5° 
(Fig.[2j , but shows quite broad wings. Also, there is the 
warp of the Galactic plane , clearly delin e ated by HI and con- 
tinuum observation s (e.g. lBurtonl fl976: Dic key fc Lockmanl 
1990; iRusseiH 12003). The importance of the warp has been 
reinforced by the latest molecular cloud surveys which re- 
veal CO clouds in the outer Galaxy mor e than 1 kpc from 
the midplane (e.g. iMizuno fc F ukui 2004). Furthermore, In- 
frared Dark Clouds (IRDCs) detected by the Midcourse 
Space Experiment (MSX), and recognis ed as potential sta r 
formation regions in their early stages l|Carev et al-l feoOO), 



A Multibeam Survey of the Galaxy for Methanol Masers 3 



1.5 



m 
£ 

0,0.5 



0) 

1 
m 



t3-0.5 

ra 

a 

-1 



■1.5 



* 



X 
X 




x 



180 



150 



120 



90 



60 30 -30 -60 

Galactic Longitude (Degrees) 



■90 



■120 



-150 -180 



Figure 1. Galactic population of 6668-MHz methanol masers within | b |sj 2°, from [Pestalozzi et alJ ll2005ll c atalogu e. Boxed regions 
highlight previous su rveys with the online v ersion ha ving blue boxes repres ent th e unbiased surveys of ICaswelll h996tJl3): green (privat e 
comm unication, see IPestalozzi et alj|2005h ; yellow jSzvmczak et al 1 l2002h ; red l lEllingsen et al.lll996l) ; purple Ipestalozzi et alj|2002h ; 
black l|Pandian et alj2007h .' Dashed line delineates the Parkes MMB survey region. Dotted line delineates the main area of the 'Piggyback' 
survey f ij5.2p . 
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Figure 2. Latitud e distribution of more than 500 methanol 
masers listed in the Pestalozzi ct al. (2005) compilation. Dotted 
line is a Gaussian fit with a FWHM of 0.52 degrees. 



have up to 50 per cent of their population outside b ^ 
0.5°. Taking these factors into account and balancing them 
against feasible observing ti me-scales, we adopted a latitude 
range of b |< 2°. For the IPestalozzi et all |2005l ) compi- 
lation 14 sources lie outside this region (four in the inner 
Galaxy, within 60 degrees of the Galactic centre, and ten in 
the outer Galaxy) . 

The velocities of the molecular clouds in which the 6668- 
MHz methanol masers reside are well traced by the emis- 



sion from the CO molecule. Thus fully sampling the ob- 
served velocity spread of CO emission is expected to ade- 
quately cover the range of velocities exhibited by 6668-MHz 
methanol masers (Section 13 . 3f) . This velocity coverage is a 
function of Galactic longitude and is shown on Fig.|3]over- 
laid on the positio n-velocity plot of the CO emission from 
iDame et al.l (|200ll ). 

The targeted la noise level of the MMB survey was 
^0.2 Jy. This is similar to the deepest previous unbiased 
surveys whic h have had la sensitiv itie s of between 0.09 Jy 
and Uy (see IPestalozzi et al.ll2005l andlPandian et alj|2007l 

!5D 



for details). From the estimates of Ivan der Waltl (|2005l ) a 3a 
threshold at <Uy should detect between about 75 and 90 
per cent of the total methanol maser population. 

P revious survey s with single dish telescopes, such as 
that of lCaswell et~afl (|l995h . have shown that often the emis- 
sion from several different methanol sources can be blended 
within the beam. Additionally, the uncertainties in the po- 
sitions derived from single dish observations at these fre- 
quencies are too large to make effective associations and to 
identify counterparts at other wavelengths. Therefore an in- 
tegral component of the MMB survey has been to determine 
the positions of the masers detected in the single dish survey 
with sub-arcsecond accuracy using interferometer observa- 
tions. 

For the Parkes detections which had not previously been 
observed at higher angular resolution, we are observing them 
with the Australia Telescope Compact Array (ATCA), and 
where possible the Multi Element Radio Linked Interferom- 
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Figure 3. Longitude versus LSR velocity (kms 1 ) plot of the greyscale of CO J = 1 — from Da me et al with grey shading 

delineating the velocity coverage of the MMB survey. 



eter Network (MERLIN); the latter, in addition to its ability 
to observe sources near declination zero and further north, 
has the advantage of even higher resolution (~15mas rms 
position uncertainty when taking into account the system- 
atic error of the phase calibration together with the individ- 
ual relative position errors). 



2.2 A Parallel Survey for Excited OH Masers 

The 1-GHz bandwidth of the MMB receiver provided the 
opportunity to survey the Galactic plane in an excited tran- 
sition of OH simultaneously with the methanol survey. Often 
seen in association with 6668-MHz methanol, the 6035-MHz 



( 2 n 3/ / 2 J = 5/2, F = 3—3) excited-state O H maser was first 



detected in the late- 1960s / earlv-1970s (e.g. l Yen et al.lll969l; 
Rvdbeck et~ai1ll970l ; IZuckerman et alJll972 ; iKnowles et all 
19731). After a numb e r of small-scale targete d observations 



i Knowles et al.lll976l; iGuilloteau et al.lll984 ). an extensive 
search by ICaswell fc Vaild (I1995T I detected 72 sources (of 
which 52 w ere new to that surv ey). Currently ~ 100 sources 
arc known (|Caswelll Il997i . l2003h . This maser's susceptibility 
to Zeeman splitting can provide additional insight into the 
magnetic fields present at sites of high-mass star formation. 



2.3 The Methanol Multibeam Receiver 

To make the MMB survey feasible required a new purpose- 
built multibeam receiver system, prompting construction 
of a new 7-beam receiver jointly by Jodrell Bank Obser- 
vatory (JBO) and the Australia Telescope National Facil- 
ity (ATNF). The receiver comprises seven dual circular po- 
larization feeds, arranged in a hexagonal pattern around 
the central feed (Fig.Q. It incorporates low noise indium- 
phosphide High Electron Mobility Transistor (HEMT) am- 
plifiers, which provide a 1-GHz bandwidth. A 500-Hz 
switched noise diode is available for continuous calibration 
in the presence of varying incoming power levels. The re- 
ceiver system is mounted at prime focus. Rotation of the 
receiver package is required to track parallalactic angle, ne- 
cessitating a complex arrangement of flexible cabling for the 
28 receiver outputs and the helium cooling lines. 




Figure 4. The footprint of the MMB receiver showing the 19.1° 
orientation of the array used to map the Galactic plane. For the 
Parkes telescope observations the footprint spans a diameter of 
approximately 15 arcmin. 



3 MULTIBEAM OBSERVATIONAL SYSTEMS 
AND ANALYSIS 

3.1 MMB Survey Techniques 

The complete MMB survey will cover the full Galactic plane, 
— 180° < / < 180° with a latitude range of ±2°. In this paper 
we focus on the southern Galactic plane portion of the sur- 
vey using the Parkes 64-m telescope; the planned northern 
hemisphere extension will adopt a similar strategy. For the 
Parkes survey region of —174° < / < 60°, the area observed 
was divided into blocks of 2° longitude by 4° latitude. Ad- 
jacent blocks have no overlap, and in the final reduction are 
combined seamlessly. The observing parameters are sum- 
marised in Table Q] 

For the Parkes MMB receiver configuration, adjacent 
beams are separated by 6.46 arcmin. The FWHM of each 
beam is 3.2 arcmin at 6668 MHz and 3.4 arcmin at 6035 
MHz. The width of the total multibeam footprint is ~15 
arcmin. For this hexagonal feed configuration, an orientation 
at an angle of 19.1° to the scan direction (the plane of the 
Galaxy, sec Fig. @ maintains equally spaced track s for the 
seven beams l|Condon. Broderick fc Sei elstad 198^). 

Scans were conducted across 2° in longitude at a rate of 
0.1° per min, dumping spectra every 5 seconds. The paral- 
lactic angle tracking was updated continuously throughout 
scans. The beams were then displaced by 1.07 arcmin in lati- 
tude before scanning back. There was then a larger displace- 
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Table 1. Table of survey parameters for the southern hemi- 
sphere observations. The positional accuracies are listed as those 
achievable with the ATCA first and where possible MERLIN sec- 
ond. 1 Noise equivalent flux density calculated from average of all 
beams and polarizations across all the observing sessions. Based 
on observations of 1934-638 and therefore applicable away from 
the Galactic plane. For the inner Galaxy there can be an addi- 
tion of several Jy from continuum emission. 2 Calculated as the 
median of the individual median noise values from each cube. 
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3.2 Targeted Follow-up 

In addition to the survey scanning observations, all detec- 
tions were later observed in pointed observations to provide 
low noise spectra. These observations tracked the source, cy- 
cling the pointing centre through each of the beams (adopt- 
ing the "MX" mode). These utilised the positions derived 
from the high-resolution observations described in Sj4] Each 
beam was on source for 1 minute, providing an effective in- 
tegration time of 7 minutes and a la noise of ~0.1 Jy. 

3.3 Multibeam data acquisition system 

Spectra were taken using the multibeam correlator and the 
wideband correlator at Parkes, integrated into a single sys- 
tem. During the first two observing runs, pending some nec- 
essary correlator modifications, data were not recorded at 
6035 MHz for the right hand circular polarization (RHCP) 
of three outer beams (2, 4, and 6). 

The observing bandwidth was 4 MHz, which corre- 
sponds to a velocity range of 180 km s -1 at 6668 MHz 
(200kms _1 at 6035 MHz ), but 150 kms" 1 was taken as 
fully usable (and we incorporated a 10 kms -1 overlap for re- 
gions with multiple centre velocities). The centre velocities 
and range for obse rvations wer e deter mined from the Galac- 
tic CO emission of lDame et all l|200ll ) (Fig.©. To fully cover 
the required velocities there was one setting for \l\ > 20°; 
two for 20° > |Z| > 6°; three for 6° > \l\ > 2°; and finally 
four for 2° > \l\ . 

The two first local oscillators were set to provide 
Doppler tracking, so that the centre of each 4-MHz band was 
kept at a constant radial velocity with respect to the kine- 
matic Local Standard of Rest (LSR) , for both the methanol 
and the OH spectral lines. The correlators were configured 
to provide 2048 frequency channels across each 4-MHz band, 
at each line frequency and each polarization. Altogether 28 
2048-channel spectra were taken every five seconds in normal 
survey scans. The channel spacing of 1950 Hz corresponds 
to a velocity spacing of 0.09 kms -1 at 6668 MHz; this pro- 
vides an effective velocity resolution of 0.11 kms -1 , which 
is sufficient to resolve the narrow maser lines (typically 
0.2— 0.3 kms -1 in width). The spectra were not smoothed. 

Fig. El shows the first light spectra obtained at Parkes 
on 22nd Jan 2006, with the source G309. 921+0.479 in the 
main beam. 

3.4 Multibeam Data Reduction 

Initial reduction of the survey data was performed with the 
program livedata and its sister program gridzilla, de- 
veloped by Mark Calabretta, for processing of multibeam 
spectral line datcQ. The capabilities of this package, and its 
first use (in the HI Parkes All Sky Survey, HIPAS S), have 
been described extensively by iBarnes et al.l (|200ll ). While 
some of our procedures are the same as those used in the 
HI survey, there are a number of major differences, as we 
describe here. In particular, we do not encounter interfer- 
ence problems, and have much less need for robust median 
statistics. 

1 http:/ /www. atnf.csiro.au/people/mcalabre/livedata. html 
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ment of 15 arcmin to the next scan pair, thus fully sampling 
a block of Galaxy of 2° in longitude by 4° in latitude with 
32 scans (16 pairs of forward and backward scans). For each 
block, a single schedule file for the Parkes Telescope Con- 
trol System (TCS) was generated to define the paths of the 
32 scans. Only one pass was needed as there was minimal 
radio frequency interference (RFI) at 6668 and 6035 MHz 
at the Parkes site. The noise equivalent flux density of the 
system temperature at both 6668 MHz and 6035 MHz was 
60 Jy. Averaged across all the observing sessions, this var- 
ied by less than 5 Jy between the different beams at each 
polarization and frequency. Across a 2° by 4° block, the 
sensitivity of a final cube generally varies by no more than 
±3.5 per cent (due to significant oversampling) , but is of 
course worse in the direction of the few continuum sources 
that are strong enough to significantly increase the system 
temperature. The pointing accuracy of the Parkes Telescope 
is typically better than 15 arcseconds. 

The noise diode is used throughout the survey observa- 
tions for amplitude calibration. However, its equivalent flux 
density value is different for each of the 28 recorded channels 
(7 beams, 2 frequencies, 2 polarizations), and to measure 
these values, and also check for pointing offsets, the contin- 
uum source 1934-638 was observed most days. This used a 
"SPOT" schedule, which scanned each beam in turn across 
the source in RA and Dec. 1934-638 was assumed to have a 
flux density of 3.9 Jy at 6668 MHz and 4.4 Jy at 6035 MHz. 
Calibration stability was further checked by observing on 
most days a known maser source using "MX" mode, in which 
a source was tracked, and the pointing centre cycled through 
each of the beams. The most commonly used source was 
G300. 969+1. 148 for both methanol and excited OH (Fig.[5] 
and [6]). From the repeatability of these measurements we es- 
timate the flux density calibration to be acccurate to a few 
percent. 

The first science run of the MMB occurred at Parkes 
between 22nd and 31st January 2006. Subsequently a further 
101 days spread over the following 22 months (split into ten 
sessions designated A to J) has yielded a total of 111 days 
of observations and completed the surveying of the Galactic 
plane between I = —174° and I = 60°. 
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Figure 5. 6668-MHz methanol source G300. 969+1. 148 as observed during observing runs 3 to 10 (for the first two runs a different source 
was observed). The median peak flux density is 4.4 Jy with maximum variation ± 10 per cent. Intrinsic source variability contributes to 
the scatter, but variability within a session was negligible. 
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Figure 6. 6035-MHz excited-state OH source G300. 969+1. 148 as observed during observing runs 3 to 10; the dashed line is the LHC 
polarization, the solid is the RHC polarization. Flux density scale refers to a single polarization. Over the period from 2006 May to 2007 
November, the LHCP had a median peak flux density of 9.0 Jy and the RHCP 6.4 Jy. The range of variation was ± 20 per cent (with 
sessions D and J showing the larg est variation from the mean) . Much of this variation is likely to be intrinsic, with the source sometimes 
twice as strong as 5 years earlier ijCaswell fc Vaildll995l) . However, similarly to the methanol, variability within a session was negligible. 
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Figure 7. First light total power spectra of G309. 921+0. 479 from 
2006 January 22. The upper panels show the 6668-MHz methanol 
spectra (right and left circular polarizations) and the lower panels 
show the 6035-MHz OH. Increasing channel number corresponds 
to increasing frequency. Details of the typical bandpass response 
can best be seen on the 6035-MHz spectra. 

A raw correlator spectrum is a composite of the maser 
spectrum, bandpass spectrum, noise and potentially, base- 
line ripples. The bandpass spectrum is the dominant compo- 
nent and was removed through the division of the correlator 
spectrum by an estimate of the bandpass (reference) spec- 
trum. 

An individual reference spectrum is required for each 
polarization of each frequency for each beam, and is ob- 
tained for each 20-minute scan using livedata. Each refer- 
ence spectrum is formed from the median of the 240 spectra 
obtained within a scan. If a maser occurs within a scan, it 
will persist for only 15 of the 240 spectra and thus have 
negligible influence on the median. The spectra were not 
smoothed since we do not have features in the spectrum 
narrower than our resolution (in contrast to the HIPASS 
survey, for which such features caused serious problems). 

Following the bandpass correction, the spectra were cal- 
ibrated for flux density using the table of values obtained 
earlier from the noise diode calibration. A fourth order poly- 
nomial was then fitted to the baseline, and the first and last 
100 edge channels filtered out. The fourth order polynomial 
provided a 15—20 per cent improvement of the systematics 
over other baseline fits, and was of a low enough order as to 
not distort any emission or absorption features that may be 
present. 

3.5 Gridding and Imaging 

The livedata processing of a 2° x 4° block of correla- 
tor data produced 32 scan files of flux density calibrated 
and bandpass corrected spectral data. Each scan consisted 
of 240 spectra for each polarization of each beam. The 
GRIDZILLA package is used to convert the individual spec- 
tra into position-position-velocity data cubes. The standard 
gridding process defines a uniform grid in the observed re- 
gion and determines how the spectra contribute to a pixel 
through applying weights which can be chosen to optimise 
various parameters, such as the signal-to-noise ratio for a 



point source, the spatial resolution, or minimising suscepti- 
bility to interference. A smoothing radius was first chosen to 
determine which spectra contribute to the pixel. The choice 
of smoothing radius affects the final gridded beamsize and 
affects the correlation between pixels and hence the image 
noise. To balance resolution and positional accuracy with 
robustness, a top-hat smoothing kernel with a diameter of 
4.4 arcmin, and hence a radius cut-off of 2.2 arcmin, was 
used. Typically, for a total intensity image, 50 spectra con- 
tribute to a pixel, 25 from each of the two polarizations. The 
pixel spacing was chosen to be 1 arcmin so as to be both a 
convenient value, and somewhat oversampled relative to the 
Nyquist requirements. 

To estimate the flux density at a given pixel, the mean of 
the contributing spectra was taken, based on the separation 
of each spectrum from the pixel. We normalise the intensity 
based on the beam response at this offset, and then down- 
weight each value by the square of the beam response. 

Initial cubes were made at the completion of each block 
(i.e. for the 32 scans a 2° x 4° image plane is constructed); 
this was adequate for initial source detection (section 13. 6p . 
but sources near the edge of a block cannot be measured pre- 
cisely without the additional information from the adjacent 
block. Since there is no overlap between adjacent observing 
blocks, in our final processing we choose to generate cubes of 
2.2° x 4°, centred at the block junctions so as to seamlessly 
blend the adjacent blocks. The increased 2.2 degree longi- 
tude span of these final cubes gave a small region of overlap 
between adjacent cubes which ensured that the analysis on 
any source could be made without any edge effects. The ve- 
locity was also cleanly limited when cubes of differing centre 
velocities were combined (Fig. [8}. 

For 6668-MHz methanol, data cubes of total intensity 
(Stokes I) were produced, combining the right and left cir- 
cular polarizations. However, for 6035-MHz OH, since the 
maser emission is typically highly polarized, additional sep- 
arate left hand circular and right hand circular polarized 
emission cubes were produced and examined. The distribu- 
tion of median noise levels across all the final cubes can be 
seen in Fig.[§] 

3.6 Source Detection 

Sources were identified in the survey cubes with an algo- 
rithm based on that of the ATNF Spectral line Analysis 
Packagfl The ASAP 1-dimension routine was run for each 
of the ~32000 spatial pixels contained within a 2.2° x 4° 
Galactic data cube. 

The algorithm operates initially on the spectrum at 
each position. The value in each spectral pixel (channel) 
is compared to a baseline estimated from a linear fit to the 
intensity as a function of frequency in a spectral window 
covering 5 per cent of the bandwidth (~ 100 channels) cen- 
tred on the channel under examination. This window was 
narrower than any baseline fluctuations, but significantly 
broader than a typical maser line. To classify as a detec- 
tion, three consecutive channels must have values more than 

4.7 times the root mean square deviation of their respective 
baseline fits (i.e. 4.7/v3 per channel). This threshold was 

2 http:/ / /www. atnf.csiro.au/computing/software/asap/ 
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Figure 8. Noise maps (after flattening the baseline) for a cube 
made with initial data (left) and final cube made from repeated 
scans (right). For the final cube the mean noise level across 24657 
pixels in the interior of the map is 0.17 Jy with a 1-sigma disper- 
sion of 0.01 Jy. The minimum value is 0.13 Jy and the maximum 
0.24 Jy and the colour scale (online version) is set to span this 
full range. The edge effect due to individual scans can be seen 
in left cube, in the area of high noise. This causes significant 
baseline fluctuation, in addition to not fully sampling and giving 
potentially erroneous positions. All are removed through the re- 
gridding process, as seen on the right, where the noisy data has 
been replaced with additional observations. 
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Figure 9. Initial distribution of noise across all the cubes for 
methanol. A few cubes were observed in poor weather, causing 
the secondary noise peak in the histogram, and these will be re- 
observed. 

found to be sufficiently high to limit the number of false 
positives (noise above the threshold), while detecting the 
majority of real sources. 

The algorithm then searched for the tails of emission by 
extending in both directions in the spectral domain until the 
difference between the channel flux density and the linear 
baseline fit changed sign. Spectral lines wholly located in the 
first 75 (lowest frequency) and last 150 (highest frequency) 
channels of each spectrum were excluded as they were likely 
to be due to the edge effects of the bandpass. Once a line was 
detected, the channels containing it were excluded and the 
baseline was recalculated. Once the end of the spectrum was 
reached the process started again, iterating until no further 
lines were found. Several lines could be detected in one pass, 



but several iterations were required to detect weak lines in 
the presence of strong lines. The noise level was determined 
at each iteration of the algorithm through calculating the 
average of the variances across the spectrum, excluding 20 
per cent of the largest absolute values as well as channels 
which met the detection criteria. 

To search for broader emission the whole process was 
repeated taking the average of each pair of channels, adopt- 
ing a smaller sample box of 50 channels and reducing the 
signal-to-noise threshold by v2. 

When a detection was made in a pixel, all neighbouring 
spatial pixels were tested in the same way, with successive 
testing of neighbouring pixels until the signal dropped below 
the threshold, or the image edge was reached. The channels 
for each pixel belonging to the new source were merged, but 
if more than one distinct spectral line component was found, 
then the spatial separation was compared. If the components 
were separated by more than one pixel (1 arcmin) the source 
was split into two (or more) appropriate positions. If sources 
were close both spatially and spectrally they were combined 
into a single source. However, there were only a few sources 
that were falsely combined and these were easily visually 
separated on inspection of the cubes. 

Once this process had iterated across every pixel in a 
cube, a source list and a series of spectra were compiled. 
The source list was verified by eye, discarding spurious de- 
tections (typically from the edge of the bandpass), and re- 
assessing the velocity range of each source. The resultant 
maser positions were accurate enough for follow-up obser- 
vations with the ATCA. The performance of the algorithm 
was tested with a cube of Gaussian noise and agreed well 
with the number of false positives expected for the chosen 
signal-to-noise threshold. The completeness of the detection 
process is discussed in detail in section [6] 

3.7 Targeted Observations Data Reduction 

The data from the targeted observations employing the 'MX' 
schedule were processed using ASAP. The bandpass reference 
for each beam was estimated using the median spectrum of 
the six 'off-source' positions, and each median reference then 
combined with the corresponding total power spectrum 'on- 
source' to determine a baseline- and gain-corrected quotient 
spectrum. For methanol the polarizations were averaged to- 
gether. The seven spectra were then combined to give a final 
spectrum with an effective 7 min integration time and a typ- 
ical rms noise of less than 0.10 Jy in the total intensity spec- 
trum. Flux densities were calibrated from the noise diode 
data and thus ultimately linked to the 'SPOT' observations 
of the continuum source 1934-638. The spectrum of a new 
detection from the survey is shown in Fig. 1101 The source 
lies in the first Galactic quadrant and its negative velocity 
shows that it is a far-side outer Galaxy source. 



4 ACCURATE POSITION MEASUREMENTS 

All southern sources detected in the Parkes MMB survey 
without a previously confirmed high resolution position were 
observed with the ATCA. Northern sources, and those near 
declination zero, were observed with MERLIN. 

The ATCA observations were taken: 2006 March/ April 
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Figure 10. An example of an outer galaxy maser found by the 
MMB survey, G027. 000-0. 041. It has kinematic Galactocentric 
and heliocentric distances of ~12kpc and ~18kpc respectively. 



and December; 2007 February, May, July and November; 
and 2008 January and August. The observations were ob- 
tained with 6-km E— W array configurations, and typically 
consisted of a loop of targets plus a nearby phase calibrator 
repeated several times within a ~10 hour period. The time 
on source was typically a few minutes, but was increased to 
more than 1 hour for weak targets. The correlator was con- 
figured with a 4- MHz bandwidth (as for the Parkes observa- 
tions) , but with 1024 frequency channels (rather than 2048) , 
so as to allow a simultaneous search for continuum emission 
at 8.64 GHz. The primary ATCA flux density calibrator 
1934-638 was observed in each session, and usually 1921-293 
as a bandpass calibrator. Processing was performed using 
the packages AIPS or miriad. Visibility data for each source 
were calibrated and then converted into a spectral cube im- 
age and cleaned. The typical rms noise is 40-60 mjy per 3.9 
kHz channel and positions are accurate to 0.1 arcsec. 

The MERLIN observations were taken between 2006 
April and 2007 November. Observations were taken in phase 
referencing mode with 5/6 cuts of 10 minutes on the source, 
at intervals of -~2 hours spread over an ~10 hour period. 
Each cut of a source was preceded and followed by ~2 min- 
utes on the phase calibrator. The observations used five of 
the six MERLIN telescopes (the MKII, Darnhall, Tabley, 
Knockin and Cambridge), with a longest baseline of 217 
km, which gives a synthesized beam size of 43 mas at 6.7 
GHz. The narrowband correlator setup was 1-MHz band- 
width and 512 channels and the wideband setup was 16- 
MHz bandwidth with 16 channels. The MERLIN data were 
initially reduced with local MERLIN software and converted 
to FITS format with preliminary b andpass scaling, the n 
all further processing was with AIPS IjDiamond et al.l l2003). 
The primary flux density and bandpass calibrator was 3C84 
(and on occasion 3C286). This was observed in both wide 
and narrowband configurations. The phase calibrator, cho- 
sen depending on the sources, was observed purely in the 
wideband setup for optimal signal-to-noise. Sources were ob- 
served at fixed frequency. Flux density and phase were cal- 
ibrated in the narrowband and wideband data respectively, 
and then the solutions applied to the source data. Inspection 
of the spectra allowed the brightest channel to be selected 



and imaged to obtain an accurate position. The typical rms 
noise is 25 mjy per channel and the positions are accurate 
to ~15mas. MERLIN observations were made of not only 
the new detections in the 20° to 60° longi tude region, but 
also t he known sources from the survey of ISzvmczak et all 
(2002) which had no published high-resolution observations. 



5 PARALLEL SURVEYS 
5.1 Magellanic Clouds Survey 

In addition to the Galactic plane survey, we have also con- 
ducted the first complete survey of the Large and Small 
Magellanic Clouds for 6668-M Hz methanol and 6035-MHz 
excited-state hydroxyl masers (|Green et al.ll200"8h . This in- 
cluded higher-sensitivity targeted searches towards known 
star-formation regions. The observations yielded the discov- 
ery of a fourth 6668-MHz methanol maser in the Large Mag- 
ellanic Cloud (LMC), found towards the star-forming re- 
gion N160a, and a second 6035-MHz excited-state hydroxyl 
maser, found towards N157a. We also re-observed the three 
previously known 6668-MHz methanol masers and the sin- 
gle 6035-MHz hydroxyl maser. All observations were initially 
made using the MMB receiver on the Parkes telescope when 
the Galactic plane was not visi ble, and accurate p ositions 
were measured with the ATCA. I Green et alj (|200ct ) present 
an analysis of the water and ground-state OH maser pop- 
ulations in addition to the excited-state OH and methanol, 
and found the LMC maser populations to be smaller than 
their Milky Way counterparts by up to a factor of ~45. We 
failed to detect emission from either transition in the Small 
Magellanic Cloud. 



5.2 Pulsar Piggyback Programme 

The 6-GHz multibeam receiver was also used for a Galac- 
tic pulsar search (Parke s project P512 c onducted by John- 
ston and colleagues; lO'Brien et al.l 1)20081 )') and we were able 
to 'piggyback' spectral line observations for the 6668-MHz 
methanol transitions on this program. The pulsar observa- 
tions were taken at two frequencies, one of which included 
the 6668-MHz methanol frequency. Piggyback observations 
occurred over eight sessions between February 2006 and Au- 
gust 2007. The requirements of the primary pulsar observa- 
tions prevented Doppler tracking and switching of the noise 
diode during the piggybacked spectral line observations. Due 
to these restrictions the methanol observations were made 
using a single fixed centre frequency (6668.5 MHz) that was 
suitable for a wide range of Galactic positions and a band- 
width of 8 MHz, with 4096 spectral channels. The piggy- 
back survey thus has the same spectral resolution as the 
main survey, but double the velocity coverage and only ob- 
served the methanol transition and not excited OH. The 
pulsar observations were made using long pointed integra- 
tions (as opposed to the scanning strategy employed for 
the MMB survey). The observations primarily targeted the 
Galactic centre, covering an area of ~51 degree 2 distributed 
over ~85 degree 2 on the sky (Fig. llip . The typical integra- 
tion times were 17 min, but a small region near the Galactic 
centre was covered by even deeper, 5-hour, integrations. 
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The different observing strategy of the piggyback obser- 
vations (compared to the main MMB survey) necessitated 
the development of a separate data reduction approach. For 
the piggyback observations the data were processed using 
ASAP. As noise diode switching was disabled for the pul- 
sar observations there were no system temperature measure- 
ments available for the piggyback data. We assumed a sys- 
tem noise equivalent flux density of 65 Jy for each beam in all 
observations, approximately 8 per cent higher than cold sky 
values (see section 3.1). We chose this higher value partly 
to be conservative in our sensitivity estimates and partly 
to account for the fact that the pulsar observations were 
more concentrated to regions of higher than average Galac- 
tic background continuum emission. The general reduction 
strategy in the case of the 17 min integration data sets was 
to form a single reference spectrum for each polarization 
and beam for each day. The reference was formed using all 
the spectra for that polarization and beam on a particular 
day and finding the median value for each channel in the 
spectrum. The quotient spectra that result from these me- 
dian references contain baseline ripples with amplitudes of 
around 1 Jy which cannot be removed even with a high order 
(> 15th) polynomial. We found that a running median fil- 
ter (of width 50 channels) produced a smoothed copy of the 
spectrum, which when subtracted from the original was very 
effective at removing the baseline ripple with minimal effect 
on real maser emission. For strong masers (peak flux den- 
sity > 10 Jy) with a velocity range comparable to the filter 
width this method results in a negative dip in the baseline 
around the maser emission. This does not effect our ability 
to detect these masers (all of which are found in the main 
MMB survey anyway) and the weak (and typically narrow 
velocity range) masers that are the primary target of the 
piggyback survey do not suffer this problem. The piggyback 
survey utilises the same line-finder algorithm employed in 
the main survey, but using a lower detection threshold of 
two consecutive channels exceeding three times the rms. 

The effective per-beam integration time for the main 
MMB survey is approximately one minute and so we would 
expect that the piggyback survey with integration times of 
17 minutes should be a factor of approximately four more 
sensitive. This is indeed the case, with the typical measured 
rms for a 17 minute integration being 42mJy. This is ap- 
proximately a factor of two more sensitive than the Arecibo 
survey of Pandian et al. (2007). Thus, the piggyback survey 
is the most sensitive large area methanol maser survey that 
has been made to date. These data provide a much deeper, 
but undersampled survey, which will be used to probe the 
low flux density end of the maser population and the vari- 
ability of these sources. Fig. Q2] shows an example spectrum 
of a weak new source from the piggyback observations; it 
demonstrates that the piggyback observations are capable 
of detecting masers with peak flux densities as low as 0.2 Jy 
and this is confirmed through numerous detections of strong 
masers in distant sidelobes. The sensitivity of the piggyback 
search is such that masers with peak flux densities in excess 
of a few hundred Jy can be weakly detected in pointings 
more than half a degree from their true location. A prelimi- 
nary analysis of approximately 25 per cent of the piggyback 
data (I ^ 10°) suggests of the order of 40 new detections, 
the majority below 0.5 Jy. 




Figure 11. Regions covered with the deep pointed integrations of 
the Piggyback survey. Shading represents actual area covered by 
beams. Solid lines outline the main region surveyed (~85 degree 2 ), 
of which 60 per cent is covered by the beams (~51 degree 2 ). 
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Figure 12. Example spectrum from the Piggyback data showing 
source G001.40+0.09. Note the low rms noise level of 0.042 Jy; this 
weak maser would not have been detected by the main survey. 

6 SURVEY SENSITIVITY & COMPLETENESS 

6.1 Efficiency of Maser Detection 

In order to test the efficiency and accuracy of the detection 
methods employed by the MMB, we constructed several test 
cubes containing simulated masers. The process involved 
creating a skymap of positions of masers together with as- 
sociated data files containing velocities and flux densities 
for each maser position. To create realistic simulated maser 
sources we took some observed 10 Jy masers and rescaled 
them to represent typical sources at a range of flux densi- 
ties. The 10 Jy masers were selected to have features across a 
realistic distribution of velocity ranges and were given ran- 
domly distributed peak velocities. The rescaled peak flux 
densities were 0.5, 0.6, 0.7, 0.8, 0.9, 1.0 and 1.1 Jy. Forty 
simulated masers were created at each flux density level, and 
split into two batches with random longitudes and latitudes 
between the ranges of 220° — 222° and —2° to +2° respec- 
tively. This range was chosen so that the simulated masers 
could be added to an empty real data cube. The simulated 
sources were then added to the bandpass calibrated data files 
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Figure 13. Results of simulation to determine detection effi- 
ciency percentage as a function of model flux density. Lines show 
the fits of the 'algorithm detections' (solid) and the final 'cata- 
logued sources' (dotted) to sigmoid functions. 



(sdfits). This was done by calculating the angular separa- 
tion and parallactic angle between the Parkes beam and the 
source position, then obtaining the attenuation factor. The 
simulated spectrum was then scaled using this factor. The 
scaled spectra were added to the data files. By adding the 
masers to the data files for a real cube, rather than simulated 
pure noise, we ensured that the effect of any residual base- 
line wobble would be taken into account. The beam mo del 
was approximated by varying the model of Ijamesl l|l987l ) to 
a dish diameter of 64 m and allowing for three feed leg shad- 
ows. This never the less had the limitation that it did not 
fully describe the sidelobe properties. 

The 14 simulated cubes (with a batch of 20 masers at a 
given flux density in each cube) were processed in the same 
manner as the survey data (section 3.7). This allowed us to 
form two sets of results: the detections made by the pro- 
gram ('algorithm detections' in Table [2] and Fig.ll3[): and 
the subset of these that were then deemed real by eye ('cat- 
alogued sources' in Tableland Fig. ll3[) . The inspection by 
eye removes any false detections by the algorithm. From the 
underlying trend we estimate the completeness to be ~80 
per cent at 0.8 Jy, approaching 100 per cent at 1 Jy. 



6.2 Completeness of Catalogue 

The recent survey with tl 

IPandian" et all (|20071 ) surveyed the region 35.2° ^ I < 53.7°, 
\b\ ^ 0.41° with a survey threshold of 0.27 Jy (more than 
95 per cent probability of detection at 3a). The MMB sur- 
vey region fully encompasses the Arecibo survey region, and 
of the 86 detections found by Pandian et al. the MMB de- 
tected 73. The strongest of the non-detections is G037.77- 
0.22 with peak flux density 0.82 Jy, and the remaining 12 
had peak flux densities (at the epoch of the Arecibo obser- 
vations) between 0.13 and 0.5 Jy. We note the Arecibo sur- 
vey had slightly narrower velocity channels, ~0.07 km s _1 , 
but was Hanning smoothed. The velocity range was —70 to 
HOkms , whilst the MMB's velocity coverage was: —35 
to 115 kms -1 for 35.2° < I < 40°; -45 to 105 kms -1 for 
40° < I < 46°; -55 to QSkms" 1 for 46° < I < 50°; and 



-65 to 85kms _1 for 50° < I < 53.7°. None of the Arecibo 
sources fall outside the MMB velocity coverage. Compar- 
isons between the two surveys are affected by the variability 
of sources from the time of the Arecibo survey to the MMB 
survey (up to three years for some sources). Some methanol 
masers exhibit large variations within less than a year (e.g . 
Caswell. Vaile fc Ellingsenll 19951: iMacLeod fc Gaviardlll996l ; 



Goedhart. Gavlard fc van der Waltll2004T r and some of the 
weaker Arecibo masers may have decreased in flux density 
below our threshold. G037.77-0.22 seems to be an example, 
since it had a recorded peak flux density of 0.82 Jy in 2005 
September, but was not detected above 0.4 Jy by the MMB 
survey in 2007 September. 

Another test of the completeness of the catalogue is ob- 
tained by taking deeper observations to much lower detec- 
tion thresholds. As described in section 5.2, we are able to 
achieve these through 'piggybacking' on pulsar observations. 
These observations will provide a better estimate of the lu- 
minosity distribution within the region and reveal what frac- 
tion of the total population we are detecting with the main 
survey. The results of this work will be fully described in a 
sub sequent paper (Ellingse n et al . in preparation). 



IPandian fc Goldsmi th (2007) in a fo l low u p paper to 



the Arecibo survey used the Ivan der Waltl (|2005l ) model dis- 
tribution to estimate that in the longitude range 40 to 50 
degrees, the fraction of total masers in the Galaxy above the 
thres hold of 0.27 Jy would be between 89 and 93 per cent. 
If the Ivan der Waltl (|2005l ) predictions are applied to MMB 
survey sensitivity we find between 79 and 86 per cent are 
above the threshold. A full analysis of the completeness by 
longitude will be discussed in future papers. 



7 SUMMARY 

The Methanol Multibeam survey, employing a new purpose- 
built 7-beam receiver, has completed observing the longitude 
range —174° < I < +60° with the Parkes radio telescope. 
More than 800 6668-MHz methanol masers have been found 
to date (including several far-side outer Galaxy masers), 
suggesting a total Galactic population in exc ess of 1000. 
We h ave also surveyed the Magellanic Clouds (|Green et al.l 
2008) and have conducted a deeper survey of a much smaller 
region in our 'piggyback' observations (Ellingsen et al. in 
preparation) . The results of the methanol survey will be re- 
leased sequentially, starting with the Galactic centre, and 
will incorporate accurate positions and spectra for each 
maser. The excited-state OH results will then be released 
and the data sets will be made available online in due course. 
Follow-up observational programmes have already been ini- 
tiated to probe the maser surroundings. 

Through incorporating new technology, multiple radio 
astronomy instruments, and international collaboration, the 
MMB survey is well on its way to producing the largest 
unbiased, homogeneous Galactic plane catalogues of 6668- 
MHz methanol and 6035-MHz hydroxyl masers to date. 
The Galactic distribution of methanol masers, when com- 
bined with distance indicators such as HI self-absorption 
and maser parallax studies, will allow for the determina- 
tion of the distribution of high-mass star formation regions, 
and thus will have implications for understanding spiral 
arm structure and Galactic dynamics. The precise source 
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Table 2. Comparison of the number of simulated sources versus those detected by the algorithm. The simulated source cube had a rms 
noise of 0.18 Jy and 40 sources were simulated at each flux density level. The first column is the input peak flux density of the simulated 
masers. The second column is the proportion of the population of simulated masers which would be expected to fall below the detection 
threshold of 3cr (0.54 Jy) based on a Gaussian distribution of flux densities about the input level. The third column is the number of 
detections made by the automated detection algorithm. The fourth column is the percentage completeness, with the associated \/N 
statistical error (again based on a Gaussian distribution). The fifth column is the number of sources deemed 'real' by visual inspection 
and the final is the completeness of this. 



Input Percentage Algorithm Completeness Catalogued Total 
Flux of sources Detections of Algorithm Sources Completeness 

(Jy) below threshold (per cent) (per cent) 



0.5 58.8 3 8±4 

0.6 36.9 11 28±8 7 18±7 

0.7 18.7 25 63±13 14 35±9 

0.8 7.4 35 88±15 34 85±15 

0.9 2.3 37 93±15 36 90±15 

1.0 sJO.l 36 90±15 32 80±14 

1.1 <0.1 40 100±16 36 90±15 



positions in the catalogues, when combined with compara- 
ble high-precision observations at sub-mm and far-IR wave- 
lengths from future instruments, will hold the key to im- 
proving our understanding of the conditions that give birth 
to high-mass stars. 
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